Gap Fluctuations in Inhomogeneous Superconductors 



O 

o 

(N 
C 



o 
o 

Or 

•4-3 



I 

C 

o 
o 



> 

On 
O 

O 



o3 



i 

c 

o 
o 



Julia S. Meyer 1 ' 2 and B. D. Simons 1 
1 Cavendish Laboratory, Madingley Road, Cambridge CB3 OHE, UK 
2 Theoretische Physik III, Ruhr-Universitat Bochum, 44780 Bochum, Germany 

(February 1, 2008) 

Spatial fluctuations of the effective pairing interaction between electrons in a superconductor 
induce variations of the order parameter which in turn lead to significant changes in the density of 
states. In addition to an overall reduction of the quasi-particle energy gap, theory suggests that 
mesoscopic fluctuations of the impurity potential induce localised tail states below the mean-field 
gap edge. Using a field theoretic approach, we elucidate the nature of the states in the 'sub-gap' 
region. Specifically, we show that these states are associated with replica symmetry broken instanton 
solutions of the mean-field equations. 

PACS numbers: 74.40. +k 74.62.Dh 74.80.-g 



I. INTRODUCTION 



The density of states (DoS) of a bulk s-wave supercon- 
ductor exhibits a quasi-particle energy gap and a singu- 
larity at the gap edge. While the gap is robust with re- 
spect to the addition of non-magnetic impurities (Ander- 
son theorem [jl|), the integrity of the gap is destroyed by 
the pair-breaking effect of time-reversal symmetry break- 
ing perturbations ||||, e.g. magnetic impurities ||, and 
parallel magnetic fields in thin films ||. (For a general 
review see, e.g., Refs. ||0]). A second, and more direct 
way of influencing the integrity of the gap is through the 
imposition of quenched spatial fluctuations of the super- 
conducting coupling constant [|]|]] . Physically, such inho- 
mogeneities can be induced by dislocations, twin or grain 
boundaries, or compositional heterogeneity as found in 
superconducting alloys jl(|. In these cases, fluctuations 
in the superconducting coupling constant are reflected in 
inhomogeneities of the order parameter which, in turn, 
induce a smearing of the quasi-particle energy gap. 

In all of the cases described above, the mechanism by 
which the quasi-particle energy gap is suppressed fol- 
lows a similar scheme and is described by the same phe- 
nomenology: at the mean-field level, each of the pertur- 
bations above lead to a suppression of the quasi-particle 
gap edge. The form of this suppression is contained 
within the Abrikosov-Gor'kov theory which describes 
the rearrangement of the ground state under the con- 
straints imposed by the self-consistency equation. While 
the physical mechanisms of gap suppression differ, the 
mean-field equations depend on a single dimensionless 
parameter characterising the strength of the external per- 
turbation (see below). Even at the mean-field level, it is 
found that if the perturbation is strong enough, the sys- 
tem is driven into a homogeneous gapless phase before 
the superconductivity is ultimately destroyed. 

After the pioneering work of Abrikosov and Gor'kov Q 
(originally formulated in the context of magnetic impu- 
rities in the dirty superconductor), it was realised that 
the integrity of the gapped phase is compromised even if 



the perturbation is weak piH13l. Optimal fluctuations 
of the random impurity potentials can conspire to cre- 
ate quasi-particle states localised on the length scale of 
the coherence length £ at energies below the mean-field 
gap, i.e. in the presence of disorder, the system fragments 
into an inhomogeneous phase in which 'droplets' of lo- 
calised sub-gap states are embedded in the superconduct- 
ing background. A similar scenario arises in proximity 
coupled systems, i.e. SN hybrid structures, where the su- 
perconductor induces a gap in the normal regio n |]P4] , |l5| . 
In the most recent investigation, it was shown ]13| that, 
close to the mean- field energy gap edge £g a p, the nature 
of the quasi-particle states (their structure and spectral 
density) are universal, depending only on the relative sep- 
aration from the edge, the dimensionality and the single 
dimensionless parameter characterising the strength of 
the perturbation. 

On this background, we consider below the influence of 
a spatially varying coupling constant g(r) on the quasi- 
particle properties of a conventional disordered s-wave 
superconductor. As mentioned above, such a program 
is not new: the same problem was investigated in an 
earlier work by Larkin and Ovchinnikov However, 
although our aims, and indeed many of our conclusions, 
are broadly similar to those of Ref. M , this investigation 
is motivated by two considerations: firstly, the develop- 
ment of a quasi-classical approach within the framework 
of the non-linear cr-model (NLcrM) to explore the nature 
of the quasi-particle states in the 'sub-gap' region serves 
as a useful prototype for future studies of related 'droplet 
phase' instabilities in other interacting theories such as 
that presented by the superconductor/insulator transi- 
tion in the disordered interacting system |p^| , |l7| . Sec- 
ondly in developing and applying the cr-model approach, 
we will find that the Lifshitz-type arguments |Q invoked 
in Ref. ||] to determine the profile of the DoS in the 
sub-gap region are flawed. Indeed, the theory developed 
below will expose a general scheme which establishes the 
universality of 'gap fluctuations' in the d-dimensional sys- 
tem (and which is in accord with the zcro-dimcnsional 
results of Ref. @). 
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With this introduction, let us formulate the model su- 
perconducting system which we will consider. Our start- 
ing point is the Gor'kov or Bogoliubov-de Gennes Hamil- 
tonian 
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-Hi 



(1) 



where the subscript 'ph' refers to the particle/hole space 
and Hq = p 2 /2m — e-p + V(r). Here, ep denotes the 
Fermi energy and V(r) represents a quenched random 
impurity potential whose distribution is characterised by 
the mean scattering time r. For simplicity, we consider 
the potential to be drawn from a Gaussian white noise 
distribution, (V(r)V(r')) v = (2TTiy T)- 1 S(r - r'), where 
vq is the average DoS of the normal system. 

The order parameter (chosen to be real) has to 
be determined self-consistcntly from the condition 
ff " 1 (r)A(r) - {^{v)^i{v))/u Q . Following Ref. §, we 
will assume that the (inverse) coupling constant g~ 1 {v) 
exhibits small fluctuations around an average value 1/g. 
As with the random impurity potential, these fluctu- 
ations of the coupling constant <5(l/<?)(r) = gi (r) are 
drawn from a Gaussian distribution with zero mean, and 
correlation 

(. 9l (r). gi (r')) 9 = ^(|r-r'|) 

(our specific choice of notation being borrowed from 
Ref. Q). Here we will impose the condition gi(r)g <C 1 so 
that the coupling constant remains positive everywhere. 
Furthermore, we will assume that the correlations are 
characterised by some correlation length r c which deter- 
mines the range of (f>. 

Qualitatively, the response of the ground state to in- 
homogeneities in the coupling constant depends sensi- 
tively on the range of the correlations. If the correlation 
length r c is much larger than the superconducting coher- 
ence length 



D 
2\A\ 



1/2 



(where D diffusion constant), then, the order parameter 
can smoothly adjust to the local value of g~ 1 (r). In this 
case A(r) ~ <?(r), and the local DoS, v{r), is fixed by the 
local value of the order parameter ||. In the opposite 
limit, one expects the faster fluctuations of the coupling 
constant to be rectified by the 'proximity effect' coupling 
of neighbouring superconducting regions. It is in this 
limit that the system becomes sensitive to quasi-classical 
phase coherence processes. Therefore, to focus our dis- 
cussion, in the following, we limit consideration to the 
quasi-classical and dirty limits, where the energy scales 
are arranged in the hierarchy 

e F > 1/t > A . 

Before turning to the formalism, let us summarise the 
main conclusions of this investigation. Firstly, as antic- 
ipated by Ref. M, it is found that inhomogeneities of 



the coupling constant are reflected in inhomogeneities of 
the superconducting order parameter. Setting A(r) = 
A + Ai(r), where A represents the homogeneous compo- 
nent of the order parameter and Ai (r) its spatial fluctu- 
ation, one finds that 

(A 1 (q)A 1 (-q)) fl = A 2 (ffi(q) fll (-q)) fl / 2 (|q|). 

Here /(|q|) represents a dimensionless function which is 
determined self-consistently (see below). 

Now, if the correlation length of the coupling constant 
is short as compared to the coherence length, r c <C £, 
one finds Q that the equation of motion for the av- 
erage quasi-particle Green function obeys a local non- 
linear equation which has the canonical form obtained 
in the theory of gapless superconductivity by Abrikosov 
and Gor'kov [Q. Specifically, in the mean-field approxi- 
mation, the BCS singularity is rounded off, and the DoS 
exhibits a reduced quasi-particle energy gap H 



2/3\3/2 



where 



£ gap = A(l - rf' A ) 



V ~ 0(0) 



£ln(£/r c ) 2 



is a dimensionless parameter characterising the strength 
of the correlations of the superconducting order param- 
eter. Note that, in the present case, the conditions 
gi(r)g < 1 (i.e. Ai(r) < A) and r c < £ imply that 
t] <C 1. Thus, as one would expect, if the coupling con- 
stant remains positive everywhere, the system remains in 
the gapped phase. 

However, the conclusions of the mean-field analysis are 
modified significantly by optimal fluctuations of the ran- 
dom impurity potential. Such fluctuations, which ap- 
pear as spatially inhomogeneous instanton configurations 
of the mean-field equation, show the gap structure to 
be fragile: mesoscopic fluctuations generate spatially lo- 
calised states at energies below the mean-field gap. Close 
to the mean- field gap edge E gap , these states are confined 
to droplets of size || 



fdrop(e) ~ £ 



A 



-1/4 



diverging as e approaches £g a p- Since r-drop £ Ap, 
each of these regions is characterised by an entire band of 
localised states confined to each droplet. To exponential 
accuracy the corresponding sub-gap DoS varies as 



v(e) ~ exp 



-a d (r 1 )u DL d - 2 ( | 
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ap 
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where a^irf) represents a known dimensionless function 
of the control parameter r/. This result, which is non- 
perturbative in the dimensionless conductance of the 
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normal system voDL d ~ 2 , differs from that obtained in 
Ref. H but, instead, mirrors the scaling obtained in the 
study of sub- gap states in superconductors with magnetic 
impurities p2| . Later we will argue that the energy scal- 
ing of the DoS is not accidental but is a universal feature 
of the sub-gap states in the superconducting system (c.f. 
Ref. 0). 

When rdrop > L, the system enters a zero-dimensional 
regime. Here the expression for the DoS (||) applies with 
d = 0. Reassuringly, this result is found to be in ac- 
cord with the exact universal result predicted for zero- 
dimensional models which exhibit a square root singular- 
ity at the level of mean-field [[lJJ . (The origin of this uni- 
versality in the present scheme was discussed in Ref. Q.) 

The paper is organised as follows: in Sec. |ll] we for- 
mulate the quantum partition function of the disordered 
superconductor in the framework of a replica field theory. 
In the dirty limit, we show that the low-energy properties 
of the bulk superconducting system is contained within 
a non-linear er-model action. Taking into account the re- 
sponse of the quasi-particles to inhomogeneities in the 
BCS coupling constant, we obtain a renormalised action 
describing gap fluctuations. Within a mean-field analysis 
we show that the quasi-particle energy gap is suppressed 
but the integrity of the gap edge is preserved. Taking 
into account instanton configurations of the action, in 
section [II we show that integrity of the gap edge is com- 
promised. An analysis of the fluctuations in the vicinity 
of the instanton configurations shows that optimal fluc- 
tuations of the impurity potential induce sub-gap states 
localised on the length scale of the superconducting co- 
herence length. Finally, in Sec. IV wc conclude on the 
universality of the results obtained here. 



II. FIELD THEORY OF THE INHOMOGENEOUS 
SUPERCONDUCTOR 



The field theory approach to the study of weakly dis- 
ordered systems |l9|-plj has been discussed and reviewed 
extensively in the literature (see, e.g., Ref. [^2)). Its ex- 
tension to the consideration of disordered superconduct- 
ing systems follows straightforwardly (2^|ij], 24| . There- 
fore, here we will only briefly summarise the main el- 
ements in the construction of the (replica) field theory 
of the disordered superconductor in the framework of 
the non-linear a-model. Using this formulation, we will 
thereafter investigate the response of the superconduct- 
ing system to inhomogeneities in the BCS coupling con- 
stant. 



A. Non-linear a-model 

The starting point of the analysis is the coherent state 
path integral for the replicated partition function [ p0|j2l] 



Z N = f f[DrDre~^ s[r K 

J a=l 

where ip a represent Grassmann fields and, defining Mat- 
subara frequencies for the Fermi system, e„ = (2n + 
1W/3, 



Howe + Sim 



with 
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dr 



drg(v) r^lVr 



To account for the symmetry properties of the system, 
it is convenient to enlarge the field space by introducing 
the four-component fields i $> aT = (tp^ , V>", ip", )/V%- 
This incorporates a particle/hole as well as time- reversal 
(or charge conjugation) space. Decoupling the quartic 
BCS interaction with the introduction of the order pa- 
rameter A(r) (chosen to be real), the total action assumes 
the canonical form 

s{* a }= f dr^2n(^n<jr~n)K 

J n 

+P °Io dT I '^W^M 

where H represents the Gor'kov Hamiltonian defined 
above ([!]). Here a^ c is a Pauli matrix in the newly intro- 
duced charge conjugation space. 

To explore the low-energy properties of the supercon- 
ducting system we follow a standard route: as with nor- 
mal disordered conductors, when subjected to an ensem- 
ble average over the random impurity distribution, the 
functional integral over the Fermionic fields '5 can be 
traded for an integral involving matrix fields Q. Physi- 
cally, the fields Q, which vary slowly on the scaling of the 
mean-free path t, describe the soft modes of density re- 
laxation — the diffusion modes. In the quasi-classical 
limit (ep 3> 1/t), the action for Q is dominated by 
the saddle-point field configuration. In the dirty limit 
(1/t ^S> A), the saddle-point equation 



Q(r) 



TW 



9 2 i i 

2^ + eF+ 27 ' |r) < 



where cr™ is a Pauli matrix in particle/hole space, admits 



Ac 



)cr5 c , with A, 



sgn(e„)(S rl 



the solution Qq - 
and [e] nm — e n 5 nm . In the limit e — > 0, the saddle point is 
not unique but spans an entire manifold parameterised 
by the unitary transformations Q — TQ$T~ l . Taking 
into account slow spatial and temporal fluctuations of 
the fields, the low-energy long-range properties of the 
weakly disordered superconducting system are described 
by a non-linear c-model action p2| , p3 17,24|] 
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S[Q,A]=v J dr J d Tff - 1 (r)A 2 (r,r) + ^ J dr tr [D(dQ) 2 - 4(ea™ ® erf + Act™)Q] , 



(3) 



where the matrix field obeys the non-linear constraint 
Q 2 = 1. The Hermitian matrices Q obey the the sym- 
metry relation Q — er™ ® crf c Q r cr™ (8) crf c reflecting the 
symmetry properties of the dyadic product \I/ £g) Note 
that the fields carry replica (a, b) as well as Matsubara 
(n, m) indices, i.e. Q = Qnm- 

This completes the formulation of the disordered su- 
perconducting system as a functional field integral in- 
volving the replicated NLerM action. Our interest here is 
in the thermodynamic DoS obtained as 



dr 

V 



r, r; ie r . 
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Making use of the analytic continuation InZ = 
limjv^o (Z N — l)/N, it is straightforward to show that 
the impurity averaged DoS can be obtained from the 
identity 



<K<0> = 



dr 

y Hm (tr [A ® ® a c 3 c QP ee }} 



Q,A 



(4) 



where (• • -)q,a = J DQ J DA - ■ ■ e~ s ^ A ^ and P e< 
projects onto the diagonal element ee. 



B. Self-consistent fluctuations of the order 
parameter 

Following Ref. [§, our strategy will be to use the mean- 
field solution of the homogeneous problem as a platform 
to develop a perturbative expansion of the self-consistent 
order parameter. Specifically, by finding the deviation 
SA(r) = Ai(r) of the order parameter from its mean 
value A to leading order in gi(r), integrating out fast 
fluctuations of Q, and averaging over random configura- 
tions of gi(r), we will obtain an effective action for the 
quasi-particle degrees of freedom of the superconduct- 
ing system. With this effective theory, we will again use 
a saddle-point analysis to explore the rearrangement of 
the ground state due to the inhomogeneous coupling con- 
stant. At the mean-field level, the solution reveals a ho- 
mogeneous renormalisation of the superconducting gap 
from its bare value. On this background, we will find 
that the hard gap predicted by the mean-field theory is 
further softened by gap fluctuations which are accom- 
modated in the effective field theory by inhomogeneous 
instanton configurations of the fields Q. 

Following this program, we begin by subjecting the 
action to a saddle-point analysis. Varying the action (||) 
with respect to Q and A, one obtains the coupled saddle- 
point equations 



Dd(QdQ) + [Q, eal« ® &. 
9 -\t)A(t) = ^tr [o?Q(t)]. 



Aaf ] = 0, 



(5a) 
(5b) 



For a homogeneous coupling constant, these equations 
admit a homogeneous solution for the order parameter 
and Q. However, for a general inhomogeneous configura- 
tion for g -1 (r), an exact solution is unavailable and an 
approximate scheme must be sought. 

Applied to the saddle-point equations (||) above, the 
Ansatz 



Qr, 



/ncr™ ® 0-3° + singer™)* 



(6) 



where 6 n = diag(0 1 l , ... , 6%) is replica diagonal, leads to 
the coupled saddle-point equations 



Dd 2 9. 



„(r) - 2 (e„ sin4(r) - A(r) cos0„(r)) = , (7a) 
1 (r)A(r) = ^sin0„(r). 



(7b) 



In the following we will drop the Matsubara indices and 
only reinstate them when necessary. These equations can 
be identified as self-consistent Usadel equations |25 26 
for the average quasi-classical Green function in the pres- 
ence of an inhomogeneous coupling constant. Specif- 
ically, the former represents the reorganisation of the 
ground state due to spatial inhomogencitics in the or- 
der parameter, while the second equation enforces the 
self-consistency condition imposed on the order parame- 
ter. 

For a homogeneous coupling constant g, the mean-field 
equations are solved by a homogeneous replica symmetric 
Ansatz with 



E 



where E 2 = e 2 + A 2 and A = (irg/p) J2 n sin #0n- In 
this expected, we simply recover the BCS solu- 

tion |Q. This result, being independent of disorder, is 
simply a manifestation of the Anderson theorem [jjj on 
the level of the effective action — in the non-interacting 
system, a weak non-magnetic impurity potential has no 
influence on the average DoS. 

To accommodate spatial fluctuations of the coupling 
constant and, with them, fluctuations of the order pa- 
rameter one should, in principle, solve the non-linear set 
of equations self-consistently. Evidently, such a program 



4 



is infeasible. Instead, following Ref. ||, taking the rela- 
tive fluctuations of the coupling constant to be small, we 
look for a perturbative expansion of the mean-field equa- 
tions. To develop the expansion, we set = 6q + #i(r) 
and, accordingly, A = A + Ai(r), where both B\ and 
Ai are of order g\. Expanding to first order in gi, one 
obtains the coupled linear equations for 6\ and Ai, 

Dd 2 6i - 2 (e cos § + A sin § )9i + 2Ai cos § = , (8a) 



E 



A x (r) = -|5]0i„(r)cos0on-ffi(r)5A. (8b) 



Transforming Eq. ( pa[ ) to the Fourier representation, one 
obtains the solution 



01 (q) =2 



Ai(q) cos §q 
L>q 2 + 2E 



which, when inserted into Eq. (8b), yields 
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FIG. 1. The function f(q), governing the dependence of 
the order parameter on variations of the coupling constant. 

Now, these fast fluctuations of the order parameter 
can have a dramatic effect on the mean-field DoS and 
its fluctuation in the vicinity of the mean-field gap edge. 
To assimilate the effect of these fluctuations it is neces- 
sary to revisit the NLerM action taking into account the 
inhomogeneous order parameter. 



Ai(q) = - 



A 



A/3 En(' 



SHI V n 



_ 2A cos 2 6 Qr 
Dq 2 +2E n 



gi (q) • 



Finally, performing the Matsubara summation, one finds 



Ax^s-A^qO/flql), 



(9) 



C. Mean- field solution 

Substituting the mean-field solution for the order pa- 
rameter, A, together with its spatial fluctuation Ai(r) 
into the NLerM action, one obtains 



S[Q] 



drtr D(dQ) 2 - 4 (£ + A^rja™ )Q 



where, normalising the wavevector q = £q by the coher- 
ence length, 



/(«) 



2q 2 



ir-{q A - l)Va In 




p + iq 4 -!) 1 / 2 

q< 1, 
q > 1. 



From this result, we obtain the response of the order 
parameter to spatial variations of the BCS coupling con- 
stant. As expected, the low Fourier components (q <C 
1 /£) of the order parameter smoothly follow spatial fluc- 
tuations of (7 -1 (r). Perhaps more surprising is the re- 
sponse of the order parameter to fast fluctuations. As 
one would expect, these fluctuations are suppressed by 
the proximity effect. However, as noted by Ref. H, the 
attenuation scales only as 1/lng 2 . 



where S = ecr™ 



Ac™. 



Now, since JdrAi(r) = 0, contributions to the gen- 
crating function arising from field configurations of Q 
which are constant or slowly varying in space are largely 
insensitive to the fluctuations. Therefore, to assess the 
influence of the spatial inhomogeneity of the order pa- 
rameter, we proceed by integrating out fast fluctuations 
of Q P7|| , where 'fast' means varying on length scales 
shorter than the coherence length. To do so, the fast and 
slow degrees are separated by expanding Q around the 
slowly varying Q, i.e. 



Q = TQT- 



Q 



-w < /2 Al 



a, 00 e w <' 2 . 



where T = exp[-IT>/2] with QIT> + IT>Q = 0. Inte- 
grating over VT>, one obtains S e s = So + S m , where 



S = 



nv 



dr tr D(dQ) 2 - 4T.Q 



and 



J 



5 



5 A 
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E 



tr 



AiCqJAxt-q'") H q ,,_ q „ [Q(q + q'), <r™] [Q(-q" - q'"), a 2 PH ] 



Here Ai is determined by Eq. (||), and n q _ q / represents the diffusion propagator, i.e. n 1 , = DcfSq^ + {£, Q(q — 
q')}. Retaining only the diagonal part n q _ q and averaging over fluctuations g±, which to a good approximation 
amounts to replacing <7i (q)<?i (— q') by its average value (<7i(q)<?i(— q'))s = 0(|q|) <5 q , q ', one obtains 



Sag 



TTVqA 

16 



dr | dr' tr [r)(r - r', e) [Q(r),af ][Q(r'), af ]] 



with 



j)(r, 



A 



(A 2 )(r)H(r). 



Upon approaching the gap edg e, n b ecomes long-ranged, 
the relevant scale being Ie — yD/E » £. Thus, the spa- 
tial dependence of 77 is governed by (A 2 ), whose range is 
determined by Eq. (0) and the correlator (gi(r)gi(r')) g ~ 
4>(\r — r'|). If <fi is short-ranged (on the scale of the co- 
herence length), we can use the approximation 



Sag — 



7W0A 
16 



drtr (7)(0,e) [Qir),al 



PH]21 



(10) 



Close to the gap, where the energy dependence of 77 is 
negligible, one has 



77 = 77 



^/^/ 2 (|q|W|qi: 



In this approximation, we will show below that the action 
recovers the mean- field equation obtained in Ref. || . 

To summarise, quenched inhomogeneities in the cou- 
pling constant induce spatial fluctuations of the order 
parameter which are accommodated by a rearrangement 
of the quasi-particles in the superconducting condensate. 
In the disordered system, this rearrangement is governed 
by the same Usadel equation that describe the proximity 
effect in hybrid SN systems. Taking into account the in- 
homogeneities in the order parameter, one obtains an ef- 
fective action for the disordered superconductor in which 
the bulk action for the non-disordered system is supple- 
mented by an additional term ( |io|) which, as we will see 
presently, leads to a suppression of the superconducting 
quasi-particle gap. 

To explore the influence of the fluctuations on the 
quasi-particle gap structure let us again vary the action 
with respect to Q. In doing so, we obtain the modified 
saddle-point equation 



Dd(QdQ)+ Q,t 



-^At; [Q, a. 



.PH 

2 



^1 



0. 



Adopting the parameterisation (|J), this saddle-point 
equation can be rewritten as 



Dd 2 9 - 2 



esinf? — Acos# + Aijsm8cos9 =0. (11 



The saddle-point equation ( pi] ) has a form which co- 
incides with the Abrikosov-Gor'kov (AG) equation ob- 
tained in the theory of gapless superconductivity [Q. 
There, the parameter 77 has to be interpreted as the spin 
scattering rate 77 = l/r s A induced by magnetic impuri- 
ties in the superconducting system. The analysis of the 
AG equation shows that, for 77 > 1 the system enters a 
gapless phase while for 77 < 1 the quasi-particle energy 
gap is suppressed but not destroyed. 

More precisely, taking the mean-field solution to be 
homogeneous in space, i.e. 9(r) = 9 AG , the saddle-point 
equation (fLl|) takes the form 



e sm ( 



Acos( 



A77 sin 6* AG cos 9 AG = 0. 



Combined with the gap equation, the solution is obtained 
self-consistently from the equation 



-r = cot ( 
A 



1 



VT + cot 2 



which coincides with Eq. (18) of Ref. [§]. Making use of 
the Eq. (|), one obtains the mean-field DoS 

u{e) = v [cos0 AG (ie„ -> e)] , 

which reveals a reduction of the energy gap according to 
i?g a p = A(l — J7 2 / 3 ) 3 / 2 . In particular, for 77 < 1, the mean- 
field theory still predicts a hard gap in the quasi-particle 
DoS, displaying square root behaviour as e — > E g&p 







7? -2/3( 1 _ ?? 2/3^-l/4 



2 e-ggap 



e<E g 
e > E, 



gap 



As expected, quenched disorder in the coupling con- 
stant is reflected in an overall suppression of the quasi- 
particle energy gap. Importantly, quasi-classical pro- 
cesses lead to a suppression of the gap that does not 
simply follow the distribution of the order parameter. 
However, as recognised by Larkin and Ovchinnikov [||, 
the square root singularity in the DoS predicted by the 
mean-field theory is untenable: optimal fluctuations as- 
sociated with the impurity potential V(r) must give rise 
to sub-gap states which cause the gap to fluctuate. How 
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are such states accommodated by the statistical field the- 
ory? 

At first sight it is tempting to seek gap fluctuations 
within the perturbative fluctuations around the symmet- 
ric mean-field saddle-point configuration 9 AG . However, 
when taken into account, it is found that the integrity 
of the gap is maintained by the analytical properties of 
the mean-field solution: perturbative fluctuations influ- 
ence only the profile of the DoS about the mean-field 
energy gap. (For a discussion of this point in the con- 
text of the hybrid SN system, see Ref. [M ). Instead, 
to explore states below the mean-field energy gap, it is 
necessary to revisit the saddle-point equation ( |ll|) and 
seek inhomogeneous 'symmetry breaking' instanton field 
configurations. 



III. INHOMOGENEOUS SADDLE POINTS: 
STRATEGY 

To develop a theory of sub-gap states in the present 
system we can draw intuition both from the analysis of 
Larkin and Ovchinnikov ||] as well as a related study 
of gap fluctuations in the superconductor with magnetic 
impurities @. In Ref. H sub- gap states were shown 
to be associated with inhomogeneous solutions of the 
Abrikosov-Gor'kov equation (O). In the framework of 
the field theory of the non-interacting system, these in- 
homogeneous instanton or bounce solutions break super- 
symmetry at the level of the action. In the present case, 
we can therefore anticipate that the relevant bounce con- 
figurations break the replica symmetry, providing an ex- 
ponential suppression of the DoS below the mean-field 
edge. At the same time, we expect to identify a zero- 
mode in the replica space which restores the global replica 
symmetry of the theory. 

In the following, for simplicity, we will first restrict at- 
tention to the quasi one-dimensional system and describe 
later how the result can be generalised to higher dimen- 
sions. 



A. Replica symmetry breaking 

In order to investigate inhomogeneous solutions, we 
first have to understand the structure of the homogeneous 
solution in more detail. The Fermionic integration con- 
tour covers the interval 9 £ [0, it}. As the mean- field DoS 
vanishes below the gap edge (i.e. v(e) — v^[cos9] = 0), 
one knows that the mean-field solution of the Abrikosov- 
Gor'kov equation satisfies the condition 9 AG (e < £gap) = 
7r/2 + i0 AG (e), where <fi AG (e) is real. 

Now instead of analysing the saddle-point equation, it 
is more convenient to study its first integral 

£ 2 (<96l) 2 + V(0) = const., 



where, switching from Matsubara to real energies, the 
effective potential is given by 



V(9) = -2i^- cos 9 + 2 sin 6+ - cos 29. 
A 2 



Along the line 9 = it/2 + i<p the potential Vr{4>) — 
— V(tt/2 + i(f>) is real with a functional dependence on 
4> shown in Fig. || for different energies. The homoge- 
neous saddle point sits at the maximum of this poten- 
tial. At e = 0, it belongs to the contour (9 AG — n/ 2) 
and the potential is symmetric around = 0. By in- 
creasing the energy, the saddle point moves away from 
the real axis along the line 9 AG = tt/2 + icf> AG until the 
energy reaches E gap . On this line the mean-field DoS 
vanishes. At energies e > -E gap the real part of 9 AG starts 
to deviate from tt/2 and the DoS becomes finite. Thus, 
even at the mean-field level the integration contour must 
be smoothly deformed to reach the saddle-point solution 
(for a discussion see, e.g., Ref. p3]). Following the be- 
haviour of the potential, one notices that one minimum 
deepens while the other becomes more and more shallow 
until merging with the maximum at e = E gap . 




-0.6 - 



-2-10 1 2 

<f 

FIG. 2. Real potential as a function of (j> along the line 
6 — 7r/2 + i(j> for e = 0, _B gap /2, -E gap . At e = -E gap the maxi- 
mum and the minimum merge. 



Now, in addition to the homogeneous Abrikosov- 
Gor'kov solution, the potential above admits for a bounce 
solution 4> AG — > </> max (> 4> AG ) — > 4>ag, where Vs(<?!> TOax ) = 
Vr(^ag)- I n principle, as is clear from Fig. ||, this is not 
the only inhomogeneous solution. However, a bounce 
solution towards negative values of ip always involves a 
larger action and its contribution can therefore be ne- 
glected. For $t[0] 7^ tt/2 the imaginary part of the poten- 
tial is finite, in general. 
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FIG. 3. Bounce solution. 

In order to obtain a finite though exponentially small 
sub-gap density of states, we need to find a solution with 
finite action. All replica symmetric solutions lead to a 
vanishing action in the limit N — > 0. Therefore, the solu- 
tion we seek necessarily involves replica symmetry break- 
ing (RSB). Leaving aside the homogeneous mean- field 
configuration, the configuration which incurs the lowest 
action is one in which the bounce inhabits only a single 
replica, say, a = 1, i.e. 



9{x) = - + i diag(>(a;), cf) AG , 



In general, an analytic solution of the highly non-linear 
saddle-point equation is not available. However, close to 
the gap edge, one can expand Vr(4>) around the homo- 
geneous solutions at E gap up to cubic order: 



V R (</>) - V R (ct> AG ) 



V6 



1/6 



A 



1/3 



where 7(e) = (E gap — e)/A and 8<p = 4> — 4>ka- In this 
limit the solution can be obtained analytically, namely 



= V6 



A 



V 2 E g 



1/6 



7 1/2 (e) cosh" 2 ( ; 



J gap / \2rdrop(^) 

Similarly, the action for the instanton assumes the form 



47tz/qA£ 



The size of the instanton, which diverges upon approach- 
ing the gap edge, is given by [0 



A 



E, 



gap 



1/12 



7 



-1/4 



(12) 



The inhomogeneous solution, thus, varies on length scales 
much longer than the coherence length justifying the sep- 
aration into slow and fast degrees of freedom in Sec. II C . 



J 



In higher dimensions one can assume that the bounce 
solution possesses radial symmetry. Nevertheless, the 
problem becomes more complicated because the saddle- 
point equation contains a gradient term P8| 



where f = |r|/£. 

However, one can still determine the parameter de- 
pendence of the action by dimensional analysis using the 
Ansatz 6<fi = a/(|r|/6). Altogether, this obtains 



Si„st. = ad(?7H^ d ~ 2 7(e) (6 - d)/4 , 



(13) 



where ad(?7) = Cd?7~ 2 / 3 (1 — ?7 2 / 3 )~( 2 + d )/ 8 and Cd a numer- 
ical constant; c\ = 2 7 7rv / 6/5. 

Although this completes our analysis of the profile and 
statistical weight of the bounce solution, since it does not 
depart from the line Jt[0] = tt/2, taken alone, it provides 
no contribution to the DoS! To understand why sub-gap 
states are associated with the bounce it is necessary to 
explore the role of fluctuations in the vicinity of the in- 
stanton. As emphasised in Ref. |Q, such a program 
turns out to be crucial in the present system. 



B. Fluctuation analysis 

As we have seen, the lowest energy bounce configura- 
tion involves a breaking of replica symmetry at the level 
of the mean-field. Taking into account fluctuations in the 
vicinity of the bounce solution, we will see below that 
there exists a zero-mode and a negative energy mode. 
The former restores global replica symmetry of the the- 
ory and, thus, ensures the integrity of the normalisation 
of the generating functional Z = 1. Furthermore, the 
negative energy mode necessitates a ir/2 rotation of the 
contour which, imparting a factor of i, renders the con- 
tribution of the instanton to the DoS non- vanishing. 

To explore the influence of the fluctuations, let us re- 
turn to the quasi one-dimensional system and introduce 
the parameterisation 



Re- w ' 2 k®o™ ®aTe w ' 2 R~ 



(14) 



where R(x) 



exp[icr™ < 



1 a^ c 8(x)/2] is the rotation from 



the metallic saddle point to the inhomogeneous saddle 
point Q S p{0). The matrices W are subject to the sym- 
metry condition in replica space Wb a — W\ b . Expanding 
the action up to second order in the generators obtains 
the following term: 



S[W] ~ J drJ2 tr [dW ab dW ba +-d<t> a d<t> b afWabPT W ba + (F a - r;cosh0 a cosh0 fc ) W ab W ba 

ab 
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-r]sinh(j)a summer™ ® tr^'Wahcr™ <g> af'Wi 



where F a = ((d</> a ) 2 + V(<f> a ) 
two types of fluctuation: 



77 cosh 20 a ) /2. There are 



(a) replica-diagonal fluctuations, and 

(b) fluctuations mixing the replicas. 

Within the replica-diagonal part, the most relevant 
contributions are due to fluctuations of the angle <fi, 
i.e. W a b = (g) a^ c ip a (x)S a b. In the N — 1 'trivial' 
replicas (9 — 9 AG ), these fluctuations are massive and, 
therefore, only lead to a weakly energy dependent pref- 
actor. More important are the fluctuations which stay 
within the symmetry broken replica a = 1 : 



S[fi 



dx I dx' if 1 (x) 



6 2 S 



5<i>{x)5<t>{x') 



fi{x') . 



Now this class of fluctuations has been studied exten- 
sively in the standard literature |29|. The operator 
5 2 S/5<fi(x)5<fi(x') has a zero- mode, <f{^ ~ d(p, due to 
translational invariance, i.e. the action is independent of 
the position of the bounce. Furthermore, as the zero- 
mode is associated with a bounce solution it has a node. 
This implies the existence of one negative energy eigen 
mode. To account for this, one has to rotate the contour 
away from the imaginary axis. This deformation of the 
contour provides a factor of i. Therefore, the result - 
which was purely imaginary ($t[6] = ir/2) before — be- 
comes real and, thus, gives a finite contribution to the 
DoS. 

Turning to the fluctuations mixing the replicas, the 
replica symmetry breaking must be accompanied by a 
zero-mode in replica space. Writing W = W~ + W + , 
where W~ {W + ) (anti-)commutes with cr™®^' , the part 
of the action involving the coupling between the replicas 
'1' and 'a 7^ 1' reads 



where 



■E 



tr 



9W? a dW± 



V[(j>{r)]-V- 



* v la * * al 



V =i (cosh20 AG + cosh2^(r) + 4cosh(</> AG ±<£(r))) . 

Although its presence is disguised by the choice of pa- 
ramcterisation, the action involving this class of fluc- 
tuations exhibits a zero-mode. Indeed, its existence is 
made manifest by specifying the parameterisation Q = 
UQ SF W. However, with this choice, the measure asso- 
ciated with the fluctuations becomes highly non-trivial. 
Nevertheless, it is useful for determining the dependence 
of the integration over zero-modes on the number of repli- 
cas, N. Close to the gap edge, i.e. at finite energies 



I 

e > 0, the structure of the saddle point within the PH- 
and CC-space is completely fixed. - The only freedom 
left are rotations in replica space, U € U(iV). Thus, 
dividing off the matrices which leave the saddle point in- 
variant, the relevant matrices U belong to the coset space 
U(iV)/(U(l)xU(iV— 1)). Integration over the zero-mode 
gives a factor which is proportional to the volume ~ N 
in the limit N — > p0| |. Or, in other words, there are N 
saddle points that contribute to the integral. Using the 
parameterisation ([l4|) does not change the N dependence 
but only the spatial structure of the zero-mode. There- 
fore, without calculating the value of the prefactor, we 
know that the result has the following form: 

v(e) ~ lim — J dx N [smh 4>{x) + {N — 1) sinh0 AG ] x 

x\ X o{x)\ 2 e- s ^- 
= dx [smhcj)(x) - sinh^ AG ] |xo(aO| 2 e~ Sinat - , 



(15) 



where | xo (f?) 1 2 describes the spatial profile of the zero- 
mode. 



C. Discussion 

This concludes our derivation of the sub-gap DoS: by 
itself the instanton or bounce configuration provides the 
leading exponential dependence of the DoS while the fluc- 
tuations render the pre-exponential factors positive defi- 
nite. More precisely, 

(1) the prefactor becomes real due to the negative en- 

ergy eigen mode and the consequential deformation 
of the contour; 

(2) the sub-gap DoS is non- vanishing only in the vicinity 

of the bounce as can be seen from Eq. (p"5"|). 

Altogether, taking the expression for the action in the 
(i-dimensional system (^3^ we obtain the expression for 
the DoS defined by Eq. (g). A further rescaling of the 
DoS allows an explicit connection between the general 
expression and its universal zero-dimensional limit: to 
this end, let us introduce the parameter 



A = ( ^A5 2 



1/3 



^4/9^2/3)1/6 



(16) 



where 6 = 1/(vqV) is the average level spacing of the 
normal system, after which the mean-field DoS in the 
vicinity of the energy gap can be cast in the form |Q 







K e > E gap ) 



E„ 



nV 



A3 



(17) 



Then, making use of Eqs. ( |16| ) and (|12|), the DoS can be 
brought to the more compact form 



Sinst. — \l — Cd 



^drop(e) 



3/2 



(18) 



where 7 g (e) = (-E gap — e )/A g and rd rop is defined in 
Eq. (12). Indeed, the factor (rd rop /L) rf can be fur- 
ther absorbed into A g by replacing the level spacing 
of the system, S, with the level spacing of the droplet, 
<Sdro P (e) = VC^ofdrop), in its definition @. Now, as 
noted in Rcf . |Q , if the mean-field DoS exhibits a square- 
root singularity of the form (|l7|) , fluctuations of the edge 
due to optimal fluctuations of the impurity potential are 
predicted to assume a universal form 



v{e) ~ exp 



E„ 



3/2 



(19) 



obtained by random matrix theory. Now, as we have 
seen, when rd TO p > L (inevitable as e — > E gap ), the sys- 
tem enters a zero-dimensional regime. In this limit, with 
d = the expression for the DoS reassuringly assumes 
the universal form (19|). 

In the present context, the mechanism by which the 
universal expression develops at the level of the action 
has been elucidated in Rcf. yjj . Specifically, in the zero- 
dimensional regime, the instanton configuration must be 
supplemented by a homogeneous replica symmetry bro- 
ken solution of stationary phase which sits at the shallow 
minimum of the potential V(0), c.f. Fig. ||. There the 
integration contour leaves the axis 9 = tt/2 + i<j), and 
the minimum represents in fact a maximum along the 
perpendicular direction . Physically, gap fluctuations 
in the zero-dimensional system correspond to sample-to- 
sample fluctuations rather than spatial variations of the 
gap. 

Although the saddle-point analysis as well as the size of 
the instanton agree with the result found in Ref. || , the 
energy dependence of the action does not, i.e. while we 
obtain the exponent a = 3/2 — d/4, the solution obtained 
in Ref. || is compatible with an exponent c*Lif = 2 — rf/4. 
As mentioned above the exponent 3/2 can be traced back 
to the zero-dimensional case, where the form of the action 
is universal 0], i.e. this exponent is a direct consequence 
of the square root behaviour of the mean- field result. 

In fact, the discrepancy of the results can be traced 
to the application of a Lifshitz-type argument to 
the present scheme. Indeed, although the problem bears 
close similarity to the Lifshitz problem of band tail states 
in semiconductors Jl^,|3l| , the correspondence is superfi- 
cial. In particular, Lifshitz tails states at the band-edge 
of a semiconductor are typically associated with wave- 
functions which vary smoothly on the scale of their ex- 
tent. As such, an estimate of the optimal character of 



the tail state distribution can be established on the level 
of the ^-field action. By contrast, the sub-gap tail states 
associated with gap fluctuations in the superconducting 
system involve a superposition of states close to the Fermi 
level, where spatial fluctuations vary rapidly oscillating 
at the scale of the Fermi wavelength — the sub-gap states 
are quasi-classical in origin. It therefore does not seem 
possible to develop a Lifshitz argument for the present 
system. As a further consequence, in contrast to the 
band-tail states, the quasi-classical nature of the sub- 
gap states in the superconductor make their properties 
insensitive to the nature of the impurity distribution. 



IV. CONCLUSIONS 



In summary, following the work of Larkin and Ovchin- 
nikov we have shown that, in a weakly disordered 
superconductor, short-scale fluctuations of the BCS cou- 
pling constant lead to a suppression of the quasi-particle 
energy gap. At the level of mean-field, the integrity of 
the gap edge is maintained. However, optimal fluctu- 
ations of the impurity potential induce a narrow band 
of states, localised at the scale of the coherence length, 
which extend below the mean-field gap edge. Within the 
framework of the statistical field theory developed here, 
these states appear as replica symmetry broken instanton 
configurations of the mean-field equations — the global 
symmetry of the theory being restored by a zero-mode in 
the replica space. To exponential accuracy, we have ob- 
tained the spectrum of gap fluctuations. The generality 
of these results has been emphasised. Specifically, in the 
c?-dimensional system, once normalised by the mean-field 
DoS at the gap edge, we have shown that the spectrum 
of tails states depends only on the dimcnsionlcss param- 
eter r\ and, in particular, is independent of the nature of 
the disorder potential. Moreover, in the zero-dimensional 
system, the spectrum of gap fluctuations is truly univer- 
sal and coincides with that obtained by Vavilov et al. 
in the study of gap fluctuations in the SN system. 

Finally, it is interesting to note that the analysis in this 
work has a number of relatives in the recent literature. As 
we have emphasised, at the level of the soft mode action, 
the theory of gap fluctuations mirrors that obtained in 
the Abrikosov-Gor'kov theory of a superconductor with 
magnetic impurities fl2| , |l3|| and, later, that encountered 
in the description of sub-gap states in the hybrid SN sys- 
tem pq| . Furthermore, various results non-perturbative 
in the (inverse) dimensionless conductance and involving 
replica (or super-)symmetry breaking have been reported 
in the literature ]3^-|35|]. Of these investigations, it is 
particularly interesting to contrast the present scheme 
with the prediction of 'anomalously localised states' in 
the weakly disordered normal conductor. 

By exploiting instanton configurations of the 
non-linear cr-model action, Khmel'nitskii and 
Muzykantskii 133] proposed that the long-time current 
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relaxation in a disordered wire was dominated by rare lo- 
calised states which coexist in a background of extended 
states (see also, Ref. These states, which are as- 

cribed to optimal fluctuations of the random potential, 
are penalised by a statistical weight which depends expo- 
nentially on the dimensionless conductance. This scaling 
mirrors that found in the present system. However, cru- 
cially, scaling in the superconducting system involves an 
energy dependence which allows the exponent to become 
small as one approaches the energy gap. 

In hindsight, it is easy to understand why optimal fluc- 
tuations can more readily induce localised states in the 
superconducting system. In the normal disordered sys- 
tem, as pointed out by Mott, hybridisation makes the co- 
existence of localised states in a background of extended 
states difficult to sustain. However, in the superconduct- 
ing system, fluctuations of the order parameter provide a 
natural mechanism by which quasi-particle states can lo- 
calise in regions where the order parameter is suppressed. 

Acknowledgements: We are grateful to Dima 
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